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Abstract 

We study the role of the 0(a)-improving clover term for the spec- 
trum of the lattice Dirac operator using cooled and thermalized SU(2) 
gauge field configurations. For cooled configurations we observe im- 
provement of the spectral properties when adding the clover term. For 
the thermalized case (12'', /3 = 2.4) without clover term we find a rather 
bad separation of physical and doubler branches making a probabilistic 
interpretation of the Atiyah-Singer index theorem on the lattice ques- 
tionable for this (3 and lattice size. Adding the clover term leads to the 
creation of additional real eigenvalues which come in pairs of opposite 
chirality thus further worsening the situation for the index theorem. 
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1. Introduction 

In the last few years the spectrum of the Dirac operator on the lattice has 
seen a lot of attention. This is partly motivated by the hope that a lattice- 
regularized gauge theory could put the idea of decomposing the fully quan- 
tized path integral into topological sectors on a conceptually sound basis. 
In the continuum gauge fields can be classified with respect to their topo- 
logical charge when they are smooth and obey certain boundary conditions. 
On the other hand the fields which carry the measure in a continuum path 
integral do not obey these conditions. The topological arguments thus can 
only be implemented on a semi-classical level. Lattice regularization might 
provide the framework to set up the topological concepts in a fully quan- 
tized theory. Up to so-called exceptional configurations lattice gauge fields 
can be assigned a topological charge and the exceptional configurations are 
expected to die out in the continuum limit. A complete decomposition of 
the path integral into topological sectors could become possible. 

In the continuum the Atiyah-Singer index theorem ||l| relates the topo- 
logical charge of classical gauge field configurations to the numbers of left- 
and right-handed zero modes of the Dirac operator. This gives rise to semi- 
classical results for fermionic observables. Again, these results can not be 
implemented in the fully quantized model since the gauge fields in the path 
integral are too rough. 

On the lattice there is no analytic result for the Atiyah-Singer index 
theorem, but it has been conjectured that on the lattice, at least for 

sufficiently smooth configurations, it is realized in a probabilistic sense. The 
question however arises if for the parameters where current simulations are 
done the gauge fields are sufficiently smooth so that the spectrum of the 
Dirac operator can be interpreted in the sense of the index theorem. 

For completely smooth configurations such as lattice approximations to 
continuum instantons there can be no doubt, that the topological interpre- 
tation of the spectrum makes sense Also for fully quantized QED2 it 
has been established |9|-[|ll|], that the index theorem governs the behavior 
of the spectrum when approaching the continuum limit. 

For thermalized lattice gauge theories in four dimensions the picture is 
less conclusive. Older results by Itoh, Iwasaki and Yoshie ||2| report a rather 
poor manifestation of the index theorem and also the distinction of physical 
and doubler modes turned out to be difficult. More recent studies [|^, p^ ] 
(compare also @],|14]-|17| for results on the spectrum for thermalized gauge 
fields) of a fermionic definition of the topological charge, which is based on 
a lattice manifestation of the index theorem, are more optimistic. However, 
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only sampled data but no detailed analysis of the index theorem for single 
thermalized configurations was presented. 

One of the motivations for this contribution is to carefully analyze the 
idea of interpreting the spectrum of the lattice Dirac operator in terms of a 
probabilistic manifestation of the Atiyah-Singer index theorem. This is done 
for a sample of thermalized SU(2) configurations on 12^ lattices at /3 = 2.4 
taken from the quenched simulation [|l8| (sample lb). The authors of [^] 
were so kind to provide us also with their cooled data obtained using the 
improved cooling method (compare also |l^]) based on the over-improved 
action given in . This allows to compare the spectra for the thermalized 
configurations with the cooled counterparts on a one by one basis. One can 
check if a fermionic definition of the topological charge based on the index 
theorem, if possible at all, agrees with the cooling procedure. Our results 
presented here indicate however, that for SU(2) quenched configurations on 
12^ lattices at /? = 2.4 most of the configurations give rise to a spectrum 
where no clear separation of physical and doubler modes can be established. 
An interpretation of the spectrum in terms of the index theorem is ques- 
tionable. 

In the last few years it was realized that improvement is a powerful if not 
essential tool when numerically analyzing lattice field theories (see e.g. [^] 
for an introductory overview). For the fermion action the 0(a)-improvement 
is governed by the clover term. Adding this extra term to the lattice Dirac 
operator certainly alters its spectral properties. 

The second motivation for this study of the eigensystem of the lattice 
Dirac operator is to see how the clover term changes the spectrum. Does it 
improve the spectral properties similar to the behavior which was observed 
in QED2 j22] using perfect actions [23]? Again we address this question by 
analyzing spectra for the thermalized, quenched sample on 12^ lattices at 
f3 = 2.4 and also for the corresponding cooled configurations. In addition 
we analyze the effect of the clover term for complete spectra using smooth 
toy configurations (constant plaquette configurations + small fiuctuations) 
on 4^ lattices. 

The article is organized as follows: In the next section we will set our 
notation and briefly discuss the framework for the Atiyah-Singer index the- 
orem on the lattice. Section 3 presents our results for the effect of the clover 
term on complete spectra using toy configurations on 4^ lattices. Section 4 
contains the study of spectra for cooled configurations on 12^ lattices. The 
index theorem and the effects of the clover term are discussed. In Section 
5 we present our results for thermalized configurations on 12^ lattices at 
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f3 = 2.4. We analyze in detail the mechanism for the proliferation of real 
eigenvalues caused by the clover term using perturbative arguments. In Sec- 
tion 6 we assess the possibility for the interpretation of the spectrum in the 
sense of the index theorem. The article closes with a discussion. 



2. Notation and basic properties of the lattice Dirac operator 

We consider the following, 0(a)-improved fermion action 

S = Y^'^{x)D{x,y)i^{y) , (1) 

x,y 

where the lattice Dirac operator D (fermion matrix) consists of three parts 

= (4 + m)I - , K = Q - Csn^C , (2) 

with m denoting the bare quark mass. For later convenience we introduced 
the abbreviation K for the non-trivial part of the fermion matrix. The 
Wilson hopping matrix Q and the clover term C are given by 



Here a^^ = i/2[7^,7i,] and F^i^ is some lattice discretization of the contin- 
uum field strength tensor. We use the standard (clover) discretization which 
can e.g. be found in |2^]. For a proper choice of the Sheikholeslami-Wohlert 
coefficient Csw in the counter-term (^), first given in |25], removes 0{a) 



cutoff effects. For SU(2) the one loop perturbative expansion gives |25, 

Cs^ = I + 0.155(1) + Oig^) , (5) 

where go is the bare coupling, related to (3 used below via /5 = 4/(7g. As in 
the case without improvement, the matrix K (or equivalently D) is neither 
hermitian nor anti-hermitian, but hermitian conjugation is implemented as 
a similarity transformation 

75K75 = . (6) 

This equation implies that the eigenvalues of K (or D) come in complex 
conjugate pairs or are real. The staggered transformation which transforms 
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Q into —Q (see e.g. [^) leaves the clover term C invariant. Thus for Cgw 7^ 
the symmetry of the spectrum of K under reflection at the imaginary axis 
is lost. 

Besides leading to symmetry of the spectrum under complex conjugation, 
(^) also has an interesting implication for the chiral properties of the 
eigenvectors of K: Only eigenvectors ip\ of K (D), where the eigenvalue A 
is real can have non-vanishing pseudoscalar matrix elements 

(V'a,75V'a) / only for A € R . (7) 

For the Dirac operator in the continuum the only eigenstates with non- 
vanishing pseudoscalar matrix elements are zero modes. Thus the eigenvec- 
tors Va of the non-hermitian matrix K (or D) which have real eigenvalues 
A should be interpreted as the "lattice zero modes" [||, ^, ^. 

Based on this interpretation of the eigenvectors with real eigenvalues as 
the lattice zero modes, one can formulate a lattice version of the Atiyah- 
Singer index theorem ^ which is expected to hold for sufficiently smooth 
gauge fields 

u[U] = R - R+ . (8) 

Here Rj^ and R- are the numbers of real eigenvalues in the physical branch 
of the spectrum with positive and negative chirality. The chirality is defined 
as the sign of the pseudoscalar matrix element (V', 75V') of the corresponding 
eigenvectors ip. i'\U\ denotes the topological charge of the lattice gauge field 
configuration. 



3. Complete spectra for test configurations on small lattices 

Relatively few is known about the role of the clover term for the spectrum 
of the lattice Dirac operator (for some results see |14, 15, |l^). To gain some 
basic insight on the role of the clover term C for all eigenvalues of K we 
computed several complete spectra of K on 4^ lattices. These computations 
were also performed to test the correct implementation of the fermion matrix 
by checking the linear and quadratic sum rules for the eigenvalues which are 
obtained by expressing the trace of K and as a functional of the gauge 
field. Both sum rules are obeyed with excellent numerical accuracy. The 
computations in this section were done using standard routines for general 
complex matrices from the LAPACK package. 

The gauge fields were constructed from constant field strength configu- 
rations where we added some fluctuations. For the details of the preparation 
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of the gauge field see [|^] . In Fig. |l] we show spectra of K for different val- 
ues of Csw The background field U has topological charge = 2 and 
is relatively smooth (in the notation of Q: s = t = 1, e = 0.1, n,. = 15). 
We show the spectra for values of the Sheikholeslami-Wohlert coefficient 
Csw = 0.0, 1.0, 1.5 and 2.5. We remark, that the gauge fields used for this 
computation are static background fields which cannot be assigned an in- 
verse squared coupling /3 and here the perturbative formula (^) can only 
give an idea about the order of magnitude of Cgw ■ 
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Figure 1: Dependence of the complete spectrum of K on Csw We plot all 
eigenvalues of K in the complex plane for a relatively smooth gauge field 
with topological charge i^[U] = 2. The coefEcient Cgw was chosen to be 
Csw = 0.0,1.0,1.5 and 2.5. There are two real eigenvalues with negative 
chirality in the physical branch (in the vicinity of 4 on the real axis), in 
accordance with (^. 
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The first plot in Fig. Q shows the spectrum of the un-improved {csw = 0) 
operator K in the complex plane. It displays the well known features of 
the spectrum for a rather smooth background gauge field with topological 
charge v = 2. The spectrum is symmetric with respect to reflection at real 
and imaginary axis. There is a bulk of complex eigenvalues, well separated 
from the real axis and several real eigenvalues clustering in small regions, 
corresponding to the physical and doubler branches. The physical branch 
(right edge of the spectrum) is well pronounced for this rather smooth gauge 
field configuration. There are no problems distinguishing the real eigenvalues 
in the physical branch of the spectrum from the doublers. The two real 
eigenvalues in the physical branch have negative chirality in accordance with 
(I)- 

When the Sheikholeslami-Wohlert coefficient Cgw is chosen equal to 1 
the spectrum changes. The symmetry under reflection at the imaginary 
axis is lost, since the clover term C is even under the staggered transfor- 
mation, while the standard hopping term Q is odd (compare the discussion 
above). However the symmetry under reflection at the real axis, which is 
due to (^) remains. The whole physical branch is shifted towards larger real 
parts (compare also the perturbative analysis in Section 5), and the complex 
eigenvalues in the physical branch seem to cluster more along a single curve. 
The discussed deformation of the spectrum becomes more pronounced when 
increasing Cgw to Csw = 1-5 and 2.5. 

It is remarkable, that for this rather smooth gauge field configuration the 
two real eigenmodes, relevant for the index theorem, remain rather unper- 
turbed when increasing Cgyj. We will later demonstrate that for thermalized 
configurations the changes of the spectrum caused by the clover term are 
much more dramatic. 

We end this section with remarking, that from the last plot it is obvious 
that the bound known at Cgw = for the eigenvalues A of -ftT: |A| < 4, or 
more general |A| < -D where D is the dimension, does not apply at Cgw > 0. 
The standard strategy [^] for proving this bound cannot be implemented 
for Csw 7^ 0, and in fact the plot shows that the bound is violated already 
for not too large Csw 



4. The physical branch of the spectrum - cooled configurations 

In this section we study in more detail the role of the improvement term 
for the physical branch of the spectrum in smooth gauge field configura- 
tions. In particular we use cooled configurations obtained from thermalized 
configurations from a quenched simulation of SU(2) gauge theory 



6 



To analyze the physical branch of the spectrum we use the Implicitly 



Restarted Arnoldi/Lanczos Method |2£] which allows to compute a few 
eigenvalues with user specified properties for large, general matrices. We 
set up the algorithm such that the 128 eigenvalues with the largest real 
parts are computed first. This gives the eigenvalues in the physical branch 
of the spectrum we are interested in. In Fig. |2| we compare the physical 
branch of the spectrum of K for Csw = and Cgw = 1 for two cooled config- 
urations from 1 18 1, one with v\U] = —1 (left hand side plot) and one with 
vlP] = 2 (right hand side plot). The value Csw = is the case without 
0(a)-improvement and Csw = 1 is the perturbative result (^ for smooth 
configurations (50 = 0). 
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Figure 2: Plots of 128 eigenvalues in the physical branch of the spectrum 
of K for cooled gauge field configurations with i'\U] = — 1 (left hand side 
plot) and i^lU] = 2 (right hand side plot). For both configurations we show 
the spectra for Cgw = (no improvement) and Cgw = 1 (perturbative value 
for smooth fields). The dashed curve is the ellipse centered at the origin 
of the complex plane with width 8 and height 4. We remark, that in the 
right-most plot the two real eigenvalues are very close to each other, such 
that the two symbols appear as one. 

The plots show, that already without improvement the complex eigen- 
values in the physical branch are concentrated in a small band along the 
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ellipse with width 8 and height 4 centered at the origin (dashed curve in 
Fig. ^). When adding the improvement term this alignment along the el- 
lipse is further enhanced. 

For the real eigenvalues the situation is different. Without improvement 
they are considerably shifted away from the limiting curve, while the clover 
term brings them very close to the ellipse, i.e. close to the origin when looking 
at the spectrum of D. Thus the real eigenvalues become exact poles of the 
improved massless propagator in smooth background fields. The number 
of real eigenvalues and their chirality is in accordance with (|8|). The value 
of the pseudoscalar matrix element is considerably improved by the clover 
term, i.e. much closer to ±1 (compare also example we give the 

values for the v = —1 configuration (left hand side plot) where we obtain 
0.8009 for the matrix element without improvement and 0.9984 at Csw = 1- 

We analyzed several more cooled configurations from sample lb of [^] 
and found the above results for cooled configurations confirmed. The num- 
ber of real eigenvalues in the physical branch and their chirality are properly 
described by ^ both for Cgw = and Cgw = 1- The clover term improves 
the alignment of the complex eigenvalues along the ellipse, shifts the real 
eigenvalues and leads to pseudoscalar matrix elements very close to ±1. 
This improvement of the spectral properties for smooth configurations when 
adding the clover term is somewhat surprising for a term which was designed 
for improving on-shell quantities. It is a little bit reminiscent of the results 
obtained for the spectrum of the perfect lattice Dirac operator in QED2 



5. Thermalized configurations 



The literature as well as the above results for the cooled con- 

figurations and the test configurations in Section 3 show, that the lattice 
index theorem (^) holds perfectly well for the spectrum in smooth gauge 
field configurations. For the case of thermalized configurations the situation 
is less clear 

In this section we carefully analyze the role of Cgw for the spectrum 
for thermalized configurations from sample lb in |15], i.e. quenched SU(2) 



gauge field configurations at /3 = 2.4 on a 12^ lattice. We compare these 
results to the results for the cooled counterparts obtained in the last section. 

Again we use the Implicitly Restarted Arnoldi/Lanczos Method to com- 
pute 200 complex eigenvalues with largest real part, i.e. eigenvalues in the 
physical branch of the spectrum. In addition we compute eigenvalues on 
the real axis using the method of counting level crossings of the hermitian 
modification of the fermion matrix. This method was first described in 101 
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and brought to further perfection and used for analyzing the real spectrum 
of the Dirac operator |, |, ||, [l|, 0. It is now also known as "overlap 



method" . 

Whenever K has a real eigenvalue r, then the operator pi —K has a 
zero eigenvalue at p = r. This implies furthermore Q, that the auxiliary 
matrix 

H{p) ^ ^^[pt- K\, (9) 

has a zero eigenvalue at p = r if and only if r is a real eigenvalue of K. 
Thus one can trace the flow of small eigenvalues p[p) of H{p) as p is varied. 
Whenever the flow crosses zero at some p = r, then this r is a real eigenvalue 
of K. From the definition (^) for the auxiliary operator H{p) it is also 
obvious that the pseudoscalar matrix element corresponding to r is given by 
the slope of the flow for the eigenvalue n{p) which crosses zero at r 0] 



(V'r,75V'r) = -1-Kp) 

dp 



(10) 



Thus the chirality assigned to r is given by the sign of the slope. The 
auxiliary operator H(p) is hermitian due to (^) and thus much simpler to 
deal with numerically. A more detailed discussion and technical remarks 
on the implementation of the level crossing algorithm can be found in the 
literature |, |, |, H |l|, ^ • 



Here we search for the real spectrum between 1.8 and 4. Real eigenval- 
ues in the vicinity of 2 are associated with the 4 doubler modes from the 
(vr, 0, 0, 0)-type corners of the Brillouin zone. In order to make sure that 
we find all real eigenvalues in the physical branch and at least some of the 
doubler s we start our search at 1.8. 

In Fig. ^ we show our results for the spectrum of K (200 complex eigen- 
values and all real eigenvalues larger than 1.8) in a thermalized background 
configuration from [|l8[ (12^ lattice, SU(2), (3 = 2.4). We remark that this is 
the thermalized configuration which corresponds to the cooled configuration 
used for the left plot in Fig. ^ (after cooling this configuration has u = — 1). 
We display the spectra for Csw = 0.0 (no improvement), Csw = 1.258 (per- 
turbative value) and Csw = 1-4 and 1.7. The last two values were chosen 
because one expects the non-perturbative value of Cgw to be larger than the 
perturbative result similar to the case of SU(3) where non-perturbative re- 
sults for Csw are known pO||- 
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Figure 3: Plots of 200 complex eigenvalues in the physical branch of the 
spectrum and all real eigenvalues larger than 1.8 for the matrix K in a 
thermalized gauge field configuration (12^ lattice, (3 = 2A). We show the 
results for Cgw = 0.0 (no improvement), Csw = 1.258 (perturbative value at 
P = 2A), for Csw = 1-4 and for Cgw = 1-7. 

An interesting feature of the plots is the strong shift of the eigenvalues 
in the physical branch of the spectrum towards larger real parts as Cg^ IS in- 
creased (compare also the perturbative arguments below) . This is equivalent 
to the known shift of the critical k when increasing Cgyj (compare e.g. the 
first article listed in [^]). We furthermore find that the eigenvalues tend to 
cluster more in a narrow band, although the improvement of the spectral 
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properties (= eigenvalues lined up along a single curve) is not as strong as 
for the cooled configurations. 

For the manifestation of the index theorem on the lattice and also other 



ideas such as the modified quenched approximation |15| the most important 
feature of the clover term is however the generation of additional real eigen- 
values. This property can be seen clearly from the four plots in Fig. ^. In 
particular we found, that the additional eigenvectors with real eigenvalues 
are generated in pairs of opposite chirality. This fact can be understood 
by analyzing how the spectral flow fi(p) for the auxiliary problem H{p) 
changes with Csw- Increasing Cgw can cause a flow line fJ.{p), which had no 
zero crossing at Cgw = 0, to develop a local maximum (or minimum) which 
eventually crosses zero from below (above) as Csw is increased further. Two 
new crossings appear, one with positive and one with negative slope. Thus, 
due to ([To|), this gives rise to two real eigenvalues with eigenvectors having 
opposite chirality. We remark that the discussed mechanism is the only way 
that additional zero crossings, i.e. real eigenvalues of K (D) can emerge, 
since the total number of zero crossings is even. This follows from the con- 
tinuity of the flow lines in p and the simple limiting behavior of H{p). For 
p — > ibcxD one has H{p) — > ^^^'^voiumexN^aior^ limiting matrices both 

have 2 x volume x Ncoior eigenvalues of +p and 2 x volume x Ncoior eigen- 
values —p. The flow lines connect the eigenvalues of the limiting cases and 
thus only an even number of zero crossings can emerge. Hence the above 
discussed mechanism for the creation of zero crossings is the only possibility 
to create new crossings, implying that the additional real eigenvalues of K 
(D) are always created in pairs of opposite chirality. 

In order to demonstrate this mechanism we isolated the flow p{p) of a 
single small eigenvalue of H{p) which exhibits the discussed phenomenon 
and present this plot in Fig. ^. We show how a particular flow line creates 
two new zero crossings as is increased from the perturbative value — 
1.258 to Csw = 1-4. We also show the flow for the intermediate values 
Csw = 1-3 and Csw = 1.35. For Csw = 1.258 the flow has no crossing of zero. 
When increasing the coefficient to Csw = 1-3 we find that a clear minimum 
has developed and the flow line is almost touching zero. At Csw = 1-35 it has 
moved further down giving rise to two crossings of zero, i.e. two additional 
real eigenvalues of K (D). At Csw = 1-4 the crossings of zero have moved 
a little and are now at approximately 2.35 (negative slope, i.e. negative 
chirality) and at 2.9 (positive chirality). 

To further study the phenomenon of proliferation of real eigenvalues 
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Figure 4: Plot of the spectral Bow fi{p) for a small eigenvalue of the auxiliary 
operator H{p) at Csw = 1-258, 1.3, 1.35 and Csw = 1-4. The flow develops 
two zero crossings as Csw is increased giving rise to two real eigenvalues of 
K (or D) which have eigenvectors with opposite chirality due to (|lO|). 



when increcLsing Cg^ we numericajlly evalua-ted. some of the mcitrix elements 
which occur in first and second order perturbation theory (in Csw) for eigen- 
values \{csw) of the operator K = Q — CgwC 



I t ^ 



^ 1 

sw) 



Kcsw) = A(0) - c^t„V'A(o)'^C'V'A(o) + cL —\7W\ + ^('^i 

(11) 

Here A(0) is the eigenvalue of K without improvement {csw = 0) and 'i/'A(o) 
is the corresponding eigenvector. Since C is a hermitian matrix it is clear 
that the first order term is real and simply shifts the eigenvalues parallel 
to the real axis. Already existing real eigenvalues remain real. This lin- 
ear term seems to be the main contribution for the shift of the spectrum 
towards larger real parts which was already discussed for Figs.||, |2| and ^. 
When analyzing the second order term for some real eigenvalue r we find 
|V'JC'0^P = iV'ICV^p. Since both, /i and /I occur in the sum for the second 
order term the imaginary parts cancel and we find that the real eigenvalue r 
is only shifted along the real axis. Again we see that existing real eigenvalues 
remain real. For complex eigenvalues A we found that the second order term 



12 



leads to a shift essentially parallel to the imaginary axis and always directed 
towards the real axis. Thus in second order complex eigenvalues can move 
closer to the real axis and eventually become real eigenvalues. 

To summarize this section we find that for thermalized gauge fields the 
clover term leads to a proliferation of real eigenvalues. The above given 
argument, based on the spectral flow for the eigenvalues of the auxiliary 
problem, shows that the eigenvalues are generated in pairs and the corre- 
sponding eigenvectors have opposite chirality. A perturbative analysis up 
to second order in Csw shows that existing real eigenvalues remain real, but 
complex eigenvalues tend to move towards the real axis. 



6. Index theorem and fermionic definition of the topological charge 

The fact that the additional real eigenvalues come in pairs with opposite chi- 
rality means that in principle they cancel each other in the lattice version 
(^) of the index theorem. However it is also obvious from Fig. ^that already 
for Csw = the separation of physical and doubler branches is not very pro- 
nounced and the additional real eigenvalues make this situation worse for 
Csw > 0. In this section we concentrate on the real spectrum and address the 
question if the separation of physical and doubler branches is large enough 
so that one can speak of a probabilistic lattice manifestation of the index 
theorem in a meaningful way. 

Fig. P shows our results for 10 thermalized configurations which for fur- 
ther reference were numbered (#1 - #10). We also quote the topological 
charges that were assigned to these configurations after cooling in [18|. The 
figure shows the real spectra for Csw = and Csw = 1.4. For both cases 
we computed all real eigenvalues larger than 1.8 in order to see at least 
a few of the real eigenvalues in the doubler branch. In the plot we also 
encode the chirality of the corresponding eigenstates. An upward pointing 
triangle means positive and a downward pointing triangle indicates negative 
chirality. 

Let's first discuss the spectra for the case without clover term. From the 
arbitrarily chosen sample we show here it is obvious, that the real eigenvalues 
are rather evenly distributed and only for configuration # 5 one can speak 
of a reasonably separated physical branch. For the other configurations it is 
unclear where to set the threshold for the physical modes. 

We conclude that for this particular setting (SU(2), 12^, (3 = 2.4) the 
naive approach to a probabilistic lattice interpretation of the Atiyah-Singer 
index theorem by simply counting the real eigenvalues in the physical branch 
seems problematic. 
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config. real spectrum at Cs„= 0.0 



real spectrum at Csw= 1-4 



#1 
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#3 
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= -1 
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#10: v = 



v= 



v= 




Figure 5: Real eigenvalues larger than 1.8 for 10 thermalized gauge field 
configurations. The value of v given on the right hand side of the plot 
was assigned after cooling |18|. An upward pointing triangle is the symbol 
for an eigenvalue with positive chirality while a downward pointing triangle 
indicates negative chirality. 



In this context it is interesting to remark, that Narayanan and Vranas 
L2| obtained a surprising result using the fermionic overlap definition of the 
topological charge and averaging over large samples of quenched gauge field 
configurations (also for SU(2), 12^, (3 = 2.4). Narayanan and Vranas take 
into account all real eigenvalues larger than 3 and discard the smaller ones 
as doublers. From our plots it is clear that even if one is willing to define 
3 to be the boundary between physical modes and doublers this definition 
disagrees in several cases 2, 6, # 7, # 8, 7^ 10) with the cooling 
results and is not justified on a one by one basis. However after sampling 
140 configurations [^], the result for the distribution of the topological 
sectors [18| from improved cooling is reproduced with good accuracy. This 
comparison of the averaged result with the analysis of single configurations 
indicates that there might be a mechanism which averages the results from 
counting the eigenvalues such that the distribution of the topological charge 
matches the outcome of the improved cooling analysis. 

As already discussed in Section 5, the clover term does not improve the 
situation for the lattice index theorem. For a larger sample this is now 
demonstrated on the right hand side of Fig. ^ where we show the 10 real 
spectra now at Sheikholeslami-Wohlert coefficient Csw = 1-4. Although the 
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additional real eigenvalues come in pairs of opposite chirality and thus in 
principle cancel in (|8|), from a technical point of view the extra real eigen- 
values are very unpleasant since the idea of dividing the eigenvalues into 
physical and doubler modes becomes even more unrealistic than at Csw = 0. 
However, also for c^^ > a phenomenon similar to the = case might 
occur, which would give the correct distribution of the topological charge 
when averaging large samples. It would be interesting to study and under- 
stand such a possibility in more detail. 

7. Discussion 

In this article we have studied the spectrum of the lattice Dirac operator 
and investigated the interplay between topological charge and spectral prop- 
erties. The effects of the 0(a)-improving clover term were analyzed. 

For relatively smooth toy configurations (constant plaquette fields -|- 
small fluctuations) on 4*^ lattices we find that the clover term leaves the 
physical branch of the spectrum rather unchanged while the bulk of the 
eigenvalues is slightly deformed. For the smooth toy conflgurations the 
lattice realization of the Atiyah-Singer index theorem is not affected by the 
clover term. 

When analyzing smooth configurations on larger lattices obtained us- 
ing the improved cooling method in [^, we even find improvement of the 
spectral properties when adding the clover term, i.e. the eigenvalues at the 
physical edge of the spectrum allign along an ellipse and the pseudoscalar 
matrix elements are close to ±1. 

For the thermalized configurations (12^,(3 = 2.4) a naive attempt to 
interpret the real spectrum in terms of a probabilistic manifestation of the 
index theorem by simply counting the real eigenvalues in the physical branch 
is problematic. We observe, that there is no reasonable separation of physical 
modes and the doubler branches. The real eigenvalues in the physical branch 
which have to be taken into account for the lattice version of the index 
theorem cannot be identified reliably for this setting. We remark however, 
that summing over large samples seems to average out the error in the 
number of real eigenvalues in the physical branch and the distribution of 
the topological charge agrees with the result from improved cooling ||l^. It 
would be interesting to understand this mechanism in more detail. 

Adding the clover term makes the separation of physical eigenmodes 
and doublers even more difficult by creating new real eigenvalues. Based 
on an argument using zero crossings of the spectral fiow of the auxiliary 
operator H{p) we showed that the additional real eigenvalues come in pairs 
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and the corresponding eigenvectors have opposite chirahty. A perturbative 
argument indicates that existing real eigenvalues are not destroyed by the 
clover term. 

Certainly there is hope, that this situation improves as one goes over to 
larger lattices and higher Such an improvement of the spectral proper- 
ties, i.e. a clearer separation of physical modes and doublers, leading to a 
simple probabilistic realization of the index theorem as one gets closer to the 
continuum limit was observed for QED2. Only further analysis can settle 
the question if a similar scenario holds for SU(N) in 4 dimensions. Also new 
approaches such as the interpolation idea [l^ or perfect actions [ 23 1 or other 



actions [31, 32] which obey the Ginsparg- Wilson relation [^] and are known 
to obey the index theorem |34] would be interesting to study numerically in 
4 dimensions. 
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